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■ Abstract. We present, for the first time, photometric and kinematical evidence, obtained with FORS2 on the VLT, for the 
existence of kinematically decoupled cores (KDCs) in two dwarf elliptical galaxies; FS76 in the NGC5044 group and FS373 

H , in the NGC3258 group. Both kinematically peculiar subcomponents rotate in the same sense as the main body of their host 

^ ' galaxy but betray their presence by a pronounced bump in the rotation velocity profiles at a radius of about 1". The KDC in 

l/") | FS76 rotates at 10 ± 3 km/s, with the host galaxy rotating at 15 ± 6 km/s; the KDC in FS373 has a rotation velocity of 6 ± 2 km/s 

while the galaxy itself rotates at 20 ± 5 km/s. FS373 has a very complex rotation velocity profile with the velocity changing 
sign at 1.5 R e . The velocity and velocity dispersion profiles of FS76 are asymmetric at larger radii. This could be caused by 
a past gravitational interaction with the giant elliptical NGC5044, which is at a projected distance of 50 kpc. We argue that 
these decoupled cores are most likely not produced by mergers in a group or cluster environment because of the prohibitively 
large relative velocities. A plausible alternative is offered by flyby interactions between a dwarf elliptical or its disky progenitor 
and a massive galaxy. The tidal forces during an interaction at the relative velocities and impact parameters typical for a group 

■ environment exert a torque on the dwarf galaxy that, according to analytical estimates, transfers enough angular momentum to 
its stellar envelope to explain the observed peculiar kinematics. 
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•'-J 1. Introduction of galaxies about 5 Gyr ago (Conselice et al. (2001) I. N-body 

rN ' simulations show that high-speed gravitational interactions 
?H ' Dwarf elliptical galaxies (dEs) are small, low-luminosity galax- . , c .. „ , , ... 

rrt • t-a j a trigger bar-formation in any small disk galaxy orbiting in 

. . . ies with diffuse, exponentially declining surface-brightness , , i A , — , j i 

_ : — -1— ,. ™ a cluster (Moore et al. (1996)1 or around a massive galaxy 

profiles (Ferguson & Binggeh (1994)). They are a gregarious 1 : . k, 1 - — ) ^nn i\t a t ■ ^ 

* i — — - J 00 in a group environment (Mayer et al. (2001)1 and strip large 

species and are found abundantly in clusters and groups of . „ . J— : — — ? r , , . 

* .... amounts of stars, gas, and dark matter from it by tidal forces, 
galaxies (although they seem to avoid the very cluster cen- T * i j ■ i «. c j- i 

Internal dynamical processes subsequently transform a disk 

ter where the tidal forces exerted by the cluster potential , . , nr. . • j i Jr .u- 

, L . galaxy into a dynamically hot sphenodal dE within a times- 
are strong enough to disrupt them (Trunllo et al. (2002))). c , . - „ c JT , . , . . . 

b b r I * i pan of about 5 Gyr. Some dEs might still contain a mem- 

According to one model for dE evolution, they are primor- ,., . , ' „ , . , , , , , 

° ory of their former state. Examples are dEs with embedded 

dial objects. Supernova explosions heat the interstellar gas . „ ,. , , , , . ir — , 

■iff o stellar disks, bars, and spiral structure (Barazza et al. (2002) 

to temperatures exceeding the escape velocity, expelling gas p .. , — , nnn , s i rrt — r ; — , hnn?\l 1 — ^ — • 

, , i T , , , * , . /' ... , bb ,, Pe Riicke et al. (2003a) Graham et al. (2003) and with siz- 

from the galaxy (Dekel& Silk (1986) Mori & Yoshii (1999)1. r., J — r Ji 7- " T~ f t - ■ 

_ . . '-r-. ; — rrr^ Jl ^tt — ^1 — a ble amounts of warm gas, suggesting recent star formation in 

This scenario explains the diffuse appearance of dEs with en- i T , , : — , , lnn ,,^ i k.. , ■ , — , ,~ n nA\i 

. r . ,.. _ , some dEs (De Ri eke et al. (2003b) Michielsen et al. (2004)). 

hanced star formation at larger radii. They are expected to ' — —. rj — 3— t— fS — 7 r t 

J . r Moreover, rotationally flattened dEs have been discovered 

form a homogeneous class and to have properties that cor- i T , . — , ,~nni\ i icr — ■ b^r ■ 1 /-i/wi\i t^, , 

. , , ., ,, . , ir T , , r (De Rijcke et al. (2001) Simien & Prugmel (2002)). The ha- 

relate tightly with mass. Alternatively, dEs could stem from ' ■ -r—. — \ 1 h ■ — ; ; -r 1 . . 

J J rassment model also offers a natural explanation for the 

late-type disk galaxies t hat entered the clusters and groups Butcher . 0emler efifect , Butcher & Ocmler ( 1 978 ) . and the 
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Fig.l. 270 sec. I-band image of FS373, a dE2,N in the 
NGC3258 group. North is up, east is left. 
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Fig. 2. 300 sec. I-band image of FS76, a dE0,N in the 
NGC5044 group. North is up, east is left. 



In this paper, we present photometric and kinematical 
evidence for the presence of kinematically decoupled cores 
(KDCs) in two dEs in a group environment: FS373 and 
FS76 (we use the galaxy identification numbers introduced 
by Fergu son & Sandage (1990)} . FS373 (Fig.Q is a nucleated 
dwarf elliptical (dE2,N) in the NGC3258 group at a distance 
of 41 Mpc (we use Ho = 70 km/s/Mpc throughout the paper). 
FS76 (Fig.E) is a dEO in the NGC5044 group, at a distance of 
36 Mpc. The pronounced bump in the rotation velocity profiles 
signals the presence of a dynamically peculiar component in 
corotation with the main body of these galaxies. Both in FS76 
and FS373, the KDC dominates the kinematics out to a radius 
of 1.5" - 2", which is well outside the nucleus or the central 
density cusp. Hence, the KDC should not be associated with 
the nucleus in the center of the host dE. It is the first time that 
evidence is found for the existence of KDCs in dwarf elliptical 
galaxies. 



Massive elliptical galaxies with KDCs both in coro- 
tation and in counterrotation with the host galaxy (e.g. 
Efst athiou et al. (1982)| |Cappellari et al. (2002)) and e llipti- 
cals with peculiar central kinematics (Fran x et al. (19 91) I have 
been known for a long time. |Bender & Su rma (1992) found 
KDCs in ellipticals to be more metalrich than the galax- 
ies' main bodies. A merger of a giant and a dwarf galaxy 
was a plausible way of producing KDCs with the observed 
properties (Balcells & Quinn (1990)1. The KDCs angular mo- 
mentum vector is set predominantly by the engulfed dwarf's 
angular momentum and need not be aligned with that of 
the giant galaxy. KDCs in SOs can equally well be ex- 
plained by the merger of two unequal-mass spiral galaxies 
(Balcells & Gonzalez (1998)1. KDCs thus provide strong ev- 
idence that mergers of (gasrich) progenitor galaxies played 
an important role in the past evolution of bright elliptical 
and lenticular galaxies, corroborating the hierarchical merging 
model for cosmological structure formation. 

In section |3 we discuss the details of the observations and 
the data-reduction process. The stellar kinematics, photome- 
try, and measurements of the strength of the near-infrared Can 
triplet (quantified by the CaT index) of these objects are pre- 
sented in section|3] followed by a study of their internal dynam- 
ics in section^] The significance of these results in the light of 
the existing theories for dE evolution is discussed in section[5] 
We summarize our conclusions in section[6] 

2. Observations and data reduction 

Within the framework of an ESO Large Program, we col- 
lected Bessel VRI-band images and deep major and minor axis 
spectra with unprecedented spatial and spectral resolution of a 
sample of 15 dEs and dwarf lenticulars (dSO), both in group 
(NGC5044, NGC5898, and NGC3258 groups) and cluster en- 
vironments (Fornax cluster). The data were taken with the 
FORS2 imaging spectrograph mounted on the VLT. The im- 
ages were bias-subtracted and flatfielded using skyflats taken 
during twilight of the same night as the science frames. The 
sky background was removed by fitting a tilted plane to regions 
of the images free of stars or other objects and subtracting it. 
The photometric zeropoints in each band were measured using 
photometric standard stars observed during the same night as 
the science frames. The images were corrected for airmass and 
interstellar extinction, using the Galactic extinction estimates 
from |Schlegel et al. (1998)| 

The spectra, with typical exposure times of 5 - 8 h per po- 
sition angle and a seeing in the range 0.3" - 1.0" FWHM, 
cover the wavelength region around the strong Can triplet 
absorption lines (~ 8600 A). All standard data reduc- 
tion procedures (bias-subtraction, flatfielding, cosmic removal, 
wavelength-calibration, sky-subtraction, flux-calibration) were 
carried out with ESO-MTDAS 1 , IRAF 2 , and our own soft- 



1 ESO-MIDAS is developed and maintained by the European 
Southern Observatory 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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ware. Fitting the dispersion relation by a cubic spline, the 
lines of the arc spectra are rectified to an accuracy of 1 - 
2 km/s FWHM. We extracted the stellar kinematical infor- 
mation by fitting a weighted mix of late G to late K gi- 
ant stars, broadened with a parameterised line-of-sight veloc- 
ity distribution (LOSVD) to the galaxy spectra. We approx- 
imated the LOSVD by a fourth-order Gauss-Hermite series 
( |Gerhard (1 993)1 1 van der Marel & Franx (1993)> (the kinemat- 
ics of the full sample will be presented in a forthcoming paper). 
The strong Can lines, which contain most of the kinematical 
information, are rather insensitive to the age and metallicity of 
an old stellar population (see Michiel sen et al. (2003)| and ref- 
erences therein), so template mismatch does not significantly 
affect the results. The spectra contain useful kinematical infor- 
mation out to 1.5-2 half-light radii (R e ). 

This is the first time a data set of dE kinematics is assem- 
bled on a par with what so far has been achieved for bright el- 
liptical galaxies. Thanks to the excellent quality of the spectra, 
both in terms of instrumental resolution and of seeing, we are 
able to spatially resolve small-scale structures in the kinematic 
profiles. 

3. Kinematics, line-strengths, and photometry 

Before discussing our observations, we first focus on the en- 
vironments of FS76 and FS373. Using ROSAT observations 
of the X-ray emitting gas in the NGC5044 group, which 
is dominated by a single central giant elliptical, NGC5044, 
|David et al. ( 1 994)1 derive a total gravitating mass of M « 
1.6 X 1O 13 /z7qM within a radius of 250hzl kpc, correspond- 
ing to M/Lb ~ 130/i5o in solar units and a galaxy veloc- 
ity dispersion <x ga i = 330 km/s. The mean systemic veloc- 
ity of the NGC5044 group is v sys = 2549 km/s, accord- 
ing to NED. The elongated NGC3258 group is dominated 
by two giant ellipticals: NGC3258 and NGC3268. ASCA ob- 
servations of the X-ray emitting gas in the NGC3258 group 
by Peder sen et al. (1997)| lead to a total mass estimate M « 
2.3 x \0 ri h~^M Q within a radius of 2A0h^ kpc, correspond- 
ing to Ml Lb ~ 240/Z50 and a galaxy velocity dispersion 
cr ga i w 400 km/s. The mean recession velocity of the group 
is v sys = 2848 km/s, according to NED. The position of FS373 
in the outskirts of the NGC3258 group is indicated in Fig. [3] 
FS76 on the other hand has a position very close to the center 
of the group (in projection), at a projected distance of 50 kpc 
west of NGC5044 (Fig.gJ. 



3.1. Kinematics 

The systemic velocity of FS373, v sys = 2415 ± 1 km/s, lies 
within 1 cr ga i of the mean recession velocity of the NGC3258 
group, confirming it as a group member. Judging from pan- 
els (b) and (d) in Fig. |5j FS373 has a complex mean velocity 
profile. The KDC's solid body rotation dominates the central 
kinematics and reaches a peak velocity of 6 ± 2 km/s around 
200 pc (or 1/8 R e ) while outside the KDC the velocity rises to 
20 ± 5 km/s at about 1 R e . The velocity changes sign again 
around 2.2 kpc * 1.5 The velocity dispersion, presented 
in panels (a) and (c) in Fig. [5] rises from 48 ± 2 km/s up to 
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Fig. 3. Position of FS373 (asterisk) inside the NGC3258 group. 
Small dots indicate the positions of the 375 group members 
listed in NED. Large dots indicate the positions of galaxies 
brighter than M B = -20. 
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Fig. 4. Position of FS76 (asterisk) inside the NGC5044 group. 
Small dots indicate the positions of the 160 group members 
listed in NED. Large dots indicate the positions of galaxies 
brighter than M B = -20. 



90 ± 12 km/s at 1 R e and remains flat beyond that radius (the 
outermost data points suggest that the dispersion may decline 
further out). Within the inner arcsecond, the velocity disper- 
sion shows a pronounced drop of about 7 km/s, suggestive of 
the presence of a dynamically cold sub-component. 

FS76 has a systemic velocity v sys = 2726 ± 2 km/s that 
places it within 1 o- ga i of the mean NGC5044 group veloc- 
ity. NGC5044 has almost the same systemic velocity: v^ 044 = 



4 



S. De Rijcke, H. Dejonghe, W. W. Zeilinger and G. K. T. Hau: dEs with KDCs 



-3.0 -2.0 



radius (kpc) 
-1.0 0.0 1.0 



2.0 3.0 



radius (kpc) 
-1.0 0.0 1.0 



in 
E 

C 

o 

'in 
i_ 

CD 
Q. 

en 



140 
120 
100 
80 
60 
40 
20 




a 

i i i I i i i i I i i i i I i i i i I i i i i I 




-15.0 -10.0 



100 



5.0 0.0 5.0 
radius (arcsec) 
radius (kpc) 
0.4 0.0 

r 



10.0 15.0 



-15.0 -10.0 



in 
E 

C 

o 

'3i 
i_ 

CD 
Q. 

cn 




5.0 0.0 5.0 
radius (arcsec) 
radius (kpc) 
0.4 0.0 

r 



10.0 15.0 



-2.0 0.0 2.0 

radius (arcsec) 




-2.0 0.0 2.0 

radius (arcsec) 



Fig. 5. Major-axis kinematics of FS373. (a) velocity dispersion, (b) mean velocity. The bottom panels zoom in on the central 
region around the KDC to show the velocity dispersion (c) and mean velocity (d) in more detail. The linear distance scale in 
kiloparsecs is indicated above the top panels (assuming Ho = 70 km/s/Mpc). The slit position on the sky is indicated in panel 
(a). FS373 shows a pronounced bump in the mean velocity profile around ~ 1" and falls to zero at ~ 2" before rising again; this 
is characteristic of a KDC. The dispersion has a central depression, within the inner arcsecond, of ~ 7 km/s, which could hint at 
the presence of a cold stellar disk. 
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Fig. 6. Major-axis kinematics of FS76. Same layout as in Fig. [5] The giant elliptical NGC5044 is at a projected distance of 50 kpc 
towards the east (i.e. positive radii are closest to NGC5044). FS76 shows a pronounced bump in the mean velocity profile around 
~ 0.8". The velocity dispersion remains approximately fiat at ~ 50 km/s within the inner 2" and rises rather sharply outside this 
radius, indicative of a cold stellar disk. 
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Table 2. R-band Sersic parameters of FS373 and FS76: extrap- 
olated central surface brightness //ro (at zero seeing), scale- 
length ro (arcsec), and shape-parameter n. 







To 




n 


FS373 


16.55 ±0.16 


0.11 ±0.02 


2.71 


±0.08 


FS76 


17.58 ±0.03 


0.50 ± 0.02 


1.85 


±0.02 



2704 + 33 km/s jNicolaci da Costa et al. (199 8)). The kinemat- 
ics of FS76, presented in Fig.|6] are equally intriguing as those 
of FS373. The major axis velocity rises rapidly to 10 ± 3 km/s 
at a radius of about 200 pc w 1/4 R e while the bulk rotation 
flattens off at 15 ± 6 km/s beyond 600 pc w 1 R e . The velocity 
dispersion is virtually constant at 50 ± 3 km/s within the inner 
2" * 0.5 R e and rises up to 75 ± 8 km/s at 1 R e . Surprisingly, 
the velocity dispersion declines beyond that radius and drops 
to 50 ± 15 km/s at 2R e . There is a hint of minor-axis rotation 
(< 5 km/s), which could be attributed to an inclined embed- 
ded disk (De Rijc ke et al. (2001)) . Outside 1 R e , the velocity 
and velocity dispersion profiles are asymmetric, suggestive of 
a past gravitational interaction, most likely with NGC5044. 

3.2. Photometry and CaT index 

We measured the surface-brightness profile, position angle, and 
ellipticity e = 1 - b/a of FS373 and FS76 as a function of 
the geometric mean of major and minor axis distance, denoted 
by a and b respectively. These were obtained using our own 
software. Basically, the code fits an ellipse through a set of po- 
sitions where a given surface brightness level is reached. The 
shape of an isophote, relative to the best fitting ellipse with 
semi-major axis a and ellipticity e, is quantified by expand- 
ing the intensity variation along this ellipse in a fourth order 
Fourier series with coefficients S4, S3, C4 and C3 : 

Ha, 6) = I (a) [1 + C 3 (a) cos(3#) + C 4 (a) cos(40)) + 

S 3 (a) sin(30)) + S 4 (a) sin(40)] . (1) 

Here, Io(a) is the average intensity of the isophote and the an- 
gle 8 is measured from the major axis. The basic photometric 
parameters of the two galaxies are presented in Table ^ The 
surface-brightness, position angle, ellipticity, and C4 profiles 
are presented in Fig. [7] 

We fitted a seeing-convolved Sersic profile to the R-band 
surface brightness profile fiR{r) = yu^o + 1.086(r/ro)^", or 
equivalently (iR(r) = //r,o + 2.5b„(r/R e ) 1,n , with ro the scale- 
length, R e the half-light radius, b„ « 0.848« - 0.142, radius 
r = yfab, measured from the central nucleus or intensity peak. 
The seeing characteristics were estimated from about 10 stars 
in each image. We simply minimized the quadratic difference 
between the observed and the seeing-convolved Sersic profile 
using a non-linear minimization routine. The formal error bars 
are approximated by the diagonal elements of the estimated co- 
variance matrix of this non-linear problem. The results of this 
fit are presented in Fig.0and Table|2] Since both galaxies have 
a central brightness peak that cannot be fitted with a Sersic pro- 
file, the inner 2" were excluded from the fit. The value n — 4 
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Fig. 8. Simulated ellipticity profile of FS373. A flat disk with 
e = 1 — b/a = 0.9 was added to a host galaxy with a flatten- 
ing varying smoothly between e = in the very center and 
e = 0.23 at 6" and beyond. This model was then convolved 
with a Gaussian to simulate FWHM=1.2" seeing conditions 
and sampled with 0.2" x 0.2" pixels in order to reproduce the 
FORS2 CCD sampling. Both the position angle and the surface 
brightness profile match the R-band characteristics of FS373. 
Without attempting to reconstruct the appearance of FS373 in 
detail, this toy model (grey full line) is able to reproduce the 
central rise of the observed ellipticity profile (black dots). 

corresponds to a de Vaucouleurs profile, typical for massive el- 
lipticals; n — 1 corresponds to the diffuse exponential profile, 
typical for the faintest dwarf ellipticals. The n values found by 
us agree with the dE classification of these galaxies. 

FS373 has a systematically positive C4 and its isophotes are 
consequently slightly disky. Moreover, the ellipticity does not 
decline all the way to the center but instead rises inwardly in- 
side the inner 2". If the true flattening were constant all the way 
to the center, convolution with a circularly symmetric Gaussian 
seeing profile would make the inner isophotes steadily rounder 
as one approaches the center. In order to understand the ob- 
served behavior of the ellipticity profile, we calculated the ap- 
pearance of a galaxy with a flattening varying smoothly be- 
tween e = in the very center and e = 2.3 at 6" and be- 
yond, to which an exponential disk with e = 9 and radial scale- 
length fiR = 0.3" is added. This model was then convolved 
with a Gaussian to simulate 1.2" FWHM seeing conditions 
and sampled with 0.2" x 0.2" pixels in order to reproduce the 
FORS2 CCD sampling. Both the position angle and the surface 
brightness profile match the R-band characteristics of FS373. 
Without attempting to reconstruct the appearance of FS373 in 
detail, this toy model (grey full line in Fig. [8} is able to repro- 
duce the central rise of the observed ellipticity profile. 

The R-I, V-R, and V-I color profiles of FS373 are pre- 
sented in Fig. (top panel). These are constructed by first fit- 
ting a cubic spline to the various surface brightness profiles as 
a function of radius r = yfab so they can all be evaluated at 
any given position. The errorbars take into account photon shot 
noise and sky-subtraction uncertainties. For each color, the im- 
age with the best seeing was convolved with a two-dimensional 
Gaussian to match the seeing of the other image (see Tabled 
for the basic photometric data). For FS373 this correction to- 
wards equal seeing was quite substantial (0.10 - 0.25 mag near 
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Table 1. Basic photometric parameters of FS76 and FS373: reddening-corrected V magnitude nty, exposure time ty (seconds), 
seeing FWHMy (arcsec), and their R and I band analogs, and the R-band effective radius R e . 





m v t v FWHMv 


m° R t R FWHM S 


m<> t, FWHM; 


Re C'/kpc) 


FS373 
FS76 


/ 720 1.5" 
15.44 180 0.6" 


14.17 360 1.2" 
14.84 220 0.9" 


13.76 270 1.1" 
14.35 300 0.8" 


7.9/1.57 
4.4/0.77 



(0 

O 




-4.0 -2.0 0.0 2.0 4.0 
radius (arcsec) 

Fig. 9. The strength of the Can triplet lines, corrected for con- 
tamination by the Pa lines, as measured by the CaT* index 
along the major axis of FS373. 
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Fig. 10. The strength of the Can triplet lines, corrected for con- 
tamination by the Pa lines, as measured by the CaT* index 
along the major axis of FS76. 



the center). The V-band image suffers from variable extinction 
by thin cirrus, affecting the surface brightness level but not the 
shape of the surface-brightness profile nor that of the color pro- 
files. Outside the central region, FS373 has slightly rising color 
profiles, an effect best seen in the V-I profile. This outward 
reddening has been observed in many dEs and is usually in- 
terpreted as a metallicity effect ( |Mori & Yo shii (1999)1. Inside 
the inner 2", all color profiles show a pronounced bump be- 
fore declining inward, although one has to take into account 



the uncertainty that comes with the large seeing correction 
that was applied. Nonetheless, this behavior correlates with 
the strength of the Can lines as measured by the CaT* index 
( |Cenarro et al. (2001)) . CaT* first rises outwards but appears to 
decline outside a radius of 2" until the index can no longer be 
measured reliably (amongst other due to residuals of the ubiq- 
uitous strong sky emission lines that make it difficult to prop- 
erly measure the continuum level). If CaT* traces metallicity 
in a composite stellar population, which is still the subject of 
debate, then one would indeed expect the KDC to be slightly 
redder than its immediate surroundings. 

Photometric parameters for FS76 are presented in Fig. 
and Tabled The isophotes of FS76 do not deviate significantly 
from ellipses (C4 w 0). Within the inner 1.2", the position an- 
gle shows a pronounced twist. Color profiles for this galaxy 
are shown in Fig. [7] The V-I color varies almost linearly with 
radius, going from V-I=1.0 mag at 0.5" to V— 1=1.1 mag at 
10". This outward reddening is also visible in the R-I profile, 
although less obvious, going from R-I=0.45 mag at 2.0" to 
R-I=0.50 mag at 10". The inner 2" appear to be slightly red- 
der again in R-I but bluer in V-R and V-I. The CaT* index 
profile, presented in Fig.E3 shows approximately the same be- 
havior as that of FS373: the strength of the Can lines rises out- 
wardly out to about 1.5" and declines further out. The higher 
CaT* value near the position where the KDC dominates the 
kinematics may either point to a slightly more metalrich (a few 
tens of a dex) or an older stellar population. Taken at face value, 
the colors and line-strengths are roughly consistent with a met- 
alpoor ([Fe/H] as -0.7) and old (T > 6 Gyr) star population. 
dVazdekis et al. (1996)||Wzdekis et al. (1997^ . 

Clearly, imaging at much higher spatial resolution, e.g. with 
HST, is required to obtain more reliable colors close to the cen- 
ter and to corroborate the presence of a disky subcomponent. 
Also, spectroscopy in a wavelength region better suited for con- 
straining the mean age and metallicity of the stellar population 
is prerequisite to fully understand the origin of these KDCs. 
Finally, it should be noted that neither galaxy shows any dust 
features. 

4. Dynamical models 

The internal dynamics of a steady-state axisymmetric stellar 
system are described by a gravitational potential \jj(m,z), that 
determines the stellar orbits, and the distribution function (DF) 
F(r, v) dr dv, which gives the number density of stars in phase 
space ((td-, z) are cylindrical coordinates). Roughly speaking, 
the DF distributes the stars over all possible orbits. According 
to the Jeans theorem, the DF can be written as a function of the 
isolating integrals of motion. An axisymmetric potential gener- 
ally allows only two such integrals; the binding energy E and 
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the z-component of the angular momentum L z . More freedom 
to distribute stars over phase space can be gained if a third inte- 
gral of motion exists. Therefore, we work in spheroidal coordi- 
nates and approximate the gravitational potential by a Stackel 
potential which allows the existence of a third integral, denoted 
by 73. Round galaxies are modeled more efficiently in a spheri- 
cal geometry, in which case the DF is conveniently taken to be 
a function of E, L z , and L, the total angular momentum. The 
L,-dependence of the DF allows the construction of rotating 
and slightly flattened stellar systems. FS373, which has a sig- 
nificantly flattened appearance on the sky, is modeled using an 
axisymmetric Stackel potential while FS76 is treated as a sys- 
tem with a spherical gravitational potential. Of course, a galaxy 
with a KDC need not be axisymmetric but in the context of 
equilibrium models such an assumption is inevitable; also, the 
isophotes do not indicate that these galaxies are significantly 
non-axisymmetric. 

A detailed account of the method we employed to construct 
the spheroidal coordinate system and Stackel potential that give 
the best fit to a given axisymmetric potential can be found 
in |Dejonghe & de Zeeuw (1988)| |Dejonghe et al. (1996)| and 
|De Bruyne*eTar*(2001)| In brief, we deproject the observed 
surface brightness distribution, derived from an I-band image, 
assuming the galaxy to be axisymmetric and viewed edge-on. 
The total mass density, including dark matter, is parameterized 
as the spatial luminosity density multiplied by a spatially vary- 
ing mass-to-light ratio 

Y( OT ,z) = a(i+B^7 2 + ( Z /<?) 2 J, (2) 

with the parameters A and B to be estimated from the data 
and q the axis ratio of the luminosity density distribution. 
The gravitational potential is obtained by decomposing the 
total mass density in spherical harmonics. Finally, we fit a 
Stackel potential to this gravitational potential. For a spher- 
ical galaxy (q = 1), the deprojection of the surface bright- 
ness profile reduces to solving an Abel's integral equation 
pSinney & Tremaine (1987) | eq. (lB-57b)). The g ravitational 
potential follows from |Binney & Tremaine (1987)1 eq. (2-22). 

For a given potential, we wish to find the DF that best re- 
produces the kinematical information. The DF is written as a 
weighted sum of basis functions F = 2,c,fi an d tne coeffi- 
cients c, are determined by minimizing the quantity 

2 v /obs, -£, <:,• obs,,, \ 2 

with obs/ an observed data point, 07 the 1-cr errorbar on that 
data point, and obs/ , the corresponding value calculated from 
the 2 th basis function, subject to the constraint that the DF 
be positive everywhere in phase space. For a spherical po- 
tential, we used the basis functions Fi(E,L) = E ai LP j and 
Fj(E,L z ) = E a 'L^ ', with integer powers, to construct the DF 
( |De Bruyne et al. (2004)) . The three-integral DF in a Stackel 
potential on the other hand was built with basis functions 
of the form Fi(E, L z , I3) = E ai L^I y J, with integer powers 
jPejonghe et al. (1996)| . In the case of FS76, the models are 
fitted directly to the observed major and minor axis spectra 



( |De Rijcke & Dejonghe (1998)) , making the best possible use 
of all the kinematical information contained in the spectra. 
The axisymmetric models for FS373 are much more complex 
and computationally time-consuming and we opted to use the 
observed surface brightness distribution and the velocity dis- 
persion and mean velocity profiles along both major and mi- 
nor axis (i.e. the central second and first order moments of 
the LOSVDs, not the Gaussian parameters) as kinematical in- 
put for the modeling code. As a further constraint, we used 
the central fourth order moment of the LOSVDs, calculated 
from the kinematic parameters up to h^. We use a Quadratic 
Programming technique to find the optimal values for the ca- 
using a library of hundreds of linearly independent basis func- 
tions, making this method essentially non-parametric. This is 
repeated for about 100 different (A, B)-pairs, covering the rele- 
vant part of parameter space. This allows us to define the range 
of models (and hence mass distributions) that are consistent 
with the data and to determine which model gives the best fit to 
the data. 

The major-axis velocity dispersion and mean velocity pro- 
files of the best models for FS373 and FS76 are compared 
with the observed kinematics in Fig. ^2 From these best mod- 
els, we estimate the mass within a 1.5 R e sphere of FS76 
at 2.9^ x 1O 9 M (at the 90% confidence level) and that of 
FS373 at 6.2^ 2 ^ x 1O 9 M . Both galaxies require the presence 
of a dark matter halo, albeit not a very massive one, in or- 
der to reproduce the observed outwardly rising velocity dis- 
persion profiles. For FS76 we find a B-band mass-to-light ra- 
tio M/Lb = 7.9tl 3 6 M Q /L B ,0 while for FS373 we find M/L B = 
7 .4^ 2 -gM / Lb,o- These models agree very well with the data 
and reproduce the central dip in the velocity dispersion, the 
bump in the velocity profile, and, in the case of FS373, the 
inward rise of the ellipticity profile. In both galaxies, this re- 
quires the presence of a fast rotating cold subcomponent in the 
models. The DFs of the best models for FS373 and FS76 are 
presented in Figs.El an dEl respectively. We plot the DF in the 
equatorial plane in turning-point space. Each orbit in this plane 
is labeled uniquely by its pericenter distance R pal and apocen- 
ter distance R- dpo if R per [ is given the same sign as L z . Circular 
orbits lie on two straight lines with R apo = ±/? pe ri- Radial orbits 
lie on the vertical line with R pei i = 0. In both galaxy mod- 
els, an excess phase-space density of stars on near-circular or- 
bits, forming the KDC, is clearly visible. Moreover, the KDC 
is obviously disjunct from the central nucleus or density cusp. 
Since the KDCs form a distinct subcomponent within their host 
galaxies, the stars that make up a KDC can be singled out of the 
DF and be studied separately (especially in the case of FS373, 
it was very clear which basis functions in the expansion of the 
DF formed the KDC). In order to roughly estimate the stellar 
mass of the KDC, we assumed a stellar mass-to-light ratio of 
M/Lb = 2 — 4M /Lb i0 , which agrees with the observed col- 
ors and line-strengths. Thus, we find Mkdc ~ 1 - 5 x 1O 7 M 
for both galaxies or a few percent at most of the total mass. 
The adopted M/Lb is both typical for a 10 Gyr old, metal-poor 
(-1 <[Fe/H]< -0.5) stellar population (which would agree 
with dEs being primordial stellar systems) and for a 5 Gyr old, 
more metal -rich (-0.5 <[Fe/H]< 0.0) stellar population (which 



8 



S. De Rijcke, H. Dejonghe, W. W. Zeilinger and G. K. T. Hau: dEs with KDCs 



radius (kpc) 
0.5 1.0 



100 




0.0 1.0 2.0 3.0 4.0 
radius (arcsec) 

radius (kpc) 
0.5 1.0 1.5 



0.0 1.0 2.0 3.0 4.0 
radius (arcsec) 

radius (kpc) 
0.2 0.4 0.6 




0.0 1.0 2.0 3.0 4.0 
radius (arcsec) 



0.0 1.0 2.0 3.0 4.0 
radius (arcsec) 



Fig. 11. Major axis velocity dispersion and mean velocity of 
FS373 (panels (a) and (c), respectively) and of FS76 (panels 
(b) and (d)), folded around the center of the galaxy. The corre- 
sponding profiles of the best fit models are plotted with a full 
line. The models for FS373 were fitted to these kinematics. On 
the other hand, the models for FS76 were fitted directly to the 
spectra and are independent from the kinematics determined 
by a Gauss-Hermite fit to the LOSVDs. Still, the model agrees 
excellently with the kinematics, as it should of course. 



would agree with dEs being harassed late-type spirals that ex- 
perienced a starburst) ( |Worthey (1994)) . 
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Fig. 12. The distribution function of FS373 in the equatorial 
plane in turning-point space. The pericenter distance, R pa i, has 
the same sign as the z-component of the angular momentum, 
L z . Circular orbits have 7? apo = ±R pa i (dashed lines); radial 
orbits are characterized by R pel -i = (dotted line). The KDC 
is visible in phase-space as an excess of stars on near-circular 
orbits, indicated by an arrow. 
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Fig. 13. The distribution function of FS76 in the equatorial 
plane in turning-point space. The KDC is visible in phase-space 
as an excess of stars on near-circular orbits. Its locus is indi- 
cated by an arrow. 



5. Discussion 

The key question is whether KDCs in dwarf elliptical galaxies 
are produced the same way as in massive ellipticals. We explore 
two possible avenues to KDC formation in dEs. The first is the 
merger hypothesis, as in giant ellipticals; the second is the ha- 
rassment scenario, which posits that gravitational interactions 
play an important role in the evolution of dEs. The analytical 
arguments given below are strictly speaking only valid for fast 
and distant encounters. An encounter between to galaxies, with 
masses Mi and M 2 , qualifies as distant if, at closest approach, 
the change in the potential energy of the pair is much smaller 
than the initial orbital kinetic energy. In a fast encounter, the 
relative velocity of the galaxies is much larger than the internal 
stellar velocities. This translates into the following constraints 
on the impact parameter b and the internal velocity dispersion 

Ota : 

2G 

b > -HM, + M 2 ); rxgai > o- int , (4) 
^1 



with V re i the relative velocity of the interacting galaxies. For 
M 2 ~ 5 x 10 9 M Q , a typical dE mass, and M\ = Mkdc « M, 
we find b > 250 - 500 pc for V re i = 0- ga i = 300 - 400 km/s. 
Also, 0-g a i > cTi nt . Hence, any non-penetrating encounter be- 
tween a dE and a much smaller dwarf galaxy classifies as a 
fast and distant encounter (even if we take the dwarf galaxy to 
be originally 10 times more massive than Mkdc, the minimum 
impact parameter would change by only 10%). In the case of a 
giant elliptical with M 2 ~ 5 x 10 n M G and Mi = M dE « M 2 , 
the condition for a fast flyby becomes b > 25 - 50 kpc, 
again rather unsensitive to Mje- In a group or cluster environ- 
ment, galaxies keep respectable distances of a few tens of kpc 
(Mo ore et al. (1996)) . With this in mind, we can discuss possi- 
ble mechanism of producing KDCs in dEs. 

5.1. The merger hypothesis 

While the merger origin of KDCs in bright ellipticals is well 
accepted, a number of facts argue against the merger hypothesis 
in the case of dEs. 



S. De Rijcke, H. Dejonghe, W. W. Zeilinger and G. K. T. Hau: dEs with KDCs 



9 



The change of the forward velocity of a galaxy with mass 
Mi induced by a fast, distant hyperbolic encounter with a 
galaxy with mass M 2 with a relative velocity V re i is given by 

2G^ 

AV^ -— r M 2 (M l + M 2 ), (5) 

( |Sparke & Gallagher (2000)1 |Binney & Tremaine (1987)) . The 
closer and the slower the encounter, the more orbital energy is 
converted into internal (stellar) kinetic energy. For an encounter 
between a typical M 2 = 5 x 1O 9 M dE and a Mj = 5x 10 7 M Q 
dwarf galaxy with a relative velocity V re i = cr ga i = 300 km/s, 
AV^j is very small (e.g. AV^| ~ 35 km/s for a collision with b = 
1 kpc). In the case of an encounter between a M 2 = 5 x 10 n M o 
elliptical and a Mi = 5 x 1O 9 M dE, on the other hand, the 
velocity change is substantial : AVy ~ V re i, even for impact 
parameters of a few tens of kiloparsecs. This suggests that a 
dE, in a group or cluster environment, has virtually no chance 
of slowing down and capturing another (smaller) dwarf galaxy, 
contrary to a more massive elliptical galaxy. Hence, once the 
galaxy group or cluster is in place, the chance of forming a 
KDC in a dE by a merger is exceedingly small. Also, it is un- 
clear how the merger scenario can explain the complex velocity 
profile of FS373, particularly the velocity changing sign around 
a radius of 12" = 2.4 kpc. 

Alternatively, the merger could have taken place before the 
group or cluster virialized, in an environment where relative 
velocities were smaller than the present values. The low galaxy 
density in such an environment argues against this idea. Also, 
it remains to be seen, e.g. using high-resolution A-body simu- 
lations, whether a KDC formed this way can survive the dE's 
falling into a group or cluster and the subsequent gravitational 
interactions with giant group or cluster members. 

5.2. The harassment scenario 

A plausible alternative is the spin-up of a dE's halo by fly-by 
encounters with other galaxies. The impulse approximation and 
the tensor virial theorem yield the following expression for the 
maximum amount of angular momentum that can be transfered 
to a galaxy with mass M\ during an encounter with a galaxy 
with mass M 2 : 

A - / - 2 W / "( 1 -?i) (6) 

with q\ the axis ratio and I\\ a component of the inertial tensor 
( |Som Sunder & Kochhar (1990)) . Using AJ ~ Mi^iAvrot to 
roughly estimate Av rot , the maximum possible change in the 
rotation velocity after re-virialization, this translates into 

with R e l the half-light radius of the galaxy with mass M\ (we 
have used I u ~ M\F? L /3). This Av 10t is plotted in Fig.^^as a 
function of impact parameter b and axis ratio q. Clearly, inter- 
actions with impact parameters small than b ~ 20 kpc can sig- 
nificantly alter the streaming velocity of these slowly rotating 
systems, especially in the case of FS373, which is apparently 
more flattened than FS76. 




10 20 30 40 50 

b (kpc) 

Fig. 14. The estimated maximum change in the rotation veloc- 
ity of a dwarf galaxy after an encounter with a giant elliptical of 
mass M 2 — 5 x 10 1 l M G as a function of impact parameter b and 
axis ratio q. Full lines: FS373 (top to bottom: q = 0.3, 0.5, 0.7), 
dotted lines: FS76 (top to bottom: q = 0.5, 0.7, 0.9). We used 
V rel = 300 km/s in the case of FS76 and V ld = 400 km/s 
for FS373. Clearly, interactions with impact parameters smaller 
than b ~ 20 kpc can signifanctly change the streaming velocity 
of these slowly rotating systems. 

According to the harassment scenario, some dEs may stem 
from late-type spirals, inflated by high-speed gravitational in- 
teractions ( |Mooreef al. (1996)| |Mayer et al. (2001)) . During 
this metamorphosis, which in the case of a late-type spiral or- 
biting a massive galaxy, e.g. in a small group, takes only 2-3 
pericenter passages, a disk and a spheroidal envelope co-exist. 
The A-body simulations presented by |Hau & Thomson (1994)| 
show that flyby interactions impart angular momentum prefer- 
entially to the outer parts of a galaxy (stars close to the galaxy's 
center have much shorter orbital periods and respond adia- 
batically to the perturbing forces during the interaction). The 
outer envelope's acquired angular momentum is most likely 
not aligned with that of the embedded disk and this results 
in a KDC-like signature in the velocity profile. Multiple in- 
teractions with different impact parameters and orbital angular 
momenta may produce the complex velocity profile observed 
in FS373, with the rotation velocity changing sign. Moreover, 
upon the first pericenter passage, the bar instability triggered 
in the progenitor spiral galaxy funnels almost all the gas to the 
center where it is consumed in a starburst. This process depletes 
the gas within roughly 2 Gyr (M ayer et al. (2001)) , effectively 
turning a gas-rich late-type galaxy into a gas-poor dE. 

6. Conclusions 

We have presented evidence for the discovery of two dwarf el- 
liptical galaxies with kinematically decoupled cores. The pres- 
ence of a cold, rotationally flattened subcomponent is sup- 
ported by the photometry, kinematics, and the dynamics of 
these galaxies. This is the first time kinematically decoupled 
cores have been detected in dwarf elliptical galaxies. 

KDCs in dwarf galaxies are not likely produced by mergers 
since the gravitational field of a dwarf elliptical galaxy is not 
strong enough for dynamical friction to sufficiently slow down 
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another dwarf galaxy at the relative velocities and impact pa- 
rameters that are typical for a group or cluster environment. In 
the field, relative velocities are much lower but the galaxy den- 
sity is prohibitively low for mergers to occur. The harassment 
scenario is able to offer an alternative explanation: the angular 
momentum, transfered to the dwarf galaxy during an encounter 
by the tidal forces, can result in the observed peculiar kinemat- 
ics. The fact that these two objects, the only ones in a sample of 
15 dEs with measured kinematics to host a KDC, are found in a 
group environment (in which mergers are unlikely but interac- 
tions are slow enough to transfer significant amounts of energy 
and angular momentum) and not in a cluster environment (in 
which mergers are unlikely and interactions are much too fast 
to transfer significant amounts of angular momentum) agrees 
with this interpretation. 
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Fig. 7. Photometry of FS373 (left column), photometry of FS76 (right column). From the bottom up: R-band surface brightness, 



Hr, overplotted with the best fitting seeing-convolved Sersic profile jr R in grey; the residual from the Sersic fit, A/i R = ji R 



/4; 

ellipticity, e = 1 - b/a; position angle, PA; and C 4 , measuring the deviations of the isophotes from a pure elliptical shape 
(these photometric quantities have all been derived from the R-band image). Vertical dotted lines indicate the seeing FWHM. 
The central nucleus, which sits in the center of the outer isophotes, has been used as zero-point for the distance scale (we use 
the geometric mean of the semi-major and semi-minor axes of each isophotal ellipse, yfab, as radius). The positivity of the C4 
parameter indicates that FS373 has disky isophotes. Top panels: R-I, V-R, and V-I color profiles as a function of major-axis 
distance a. For each color index, the image with the best seeing was convolved with a Gaussian to match the seeing of the other 
image (vertical dotted lines indicate where the seeing correction becomes important, with the V-R and V-I colors having a worse 
seeing (grey line) than the R-I color (black line) in the case of FS373). In the case of FS373, the V-band image suffers from 
variable extinction by thin cirrus, affecting the V-band surface brightness level but not the shape of the surface-brightness profile 
nor that of the color profiles. The uncorrected color profiles are drawn in full lines. The seeing correction is quite substantial, 
especially in the profiles involving the V-band image (0.2-0.3 mag). For FS76, the seeing corrections are small; 0.05 at most. 
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